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Infrared Features of the Lattice Landau Gauge QCD 
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The infrared properties of lattice Landau gauge QCD of S!7(3) are studied by measuring gluon propagator, 
ghost propagator, QCD running coupling and Kugo-Ojima parameter of (3 = 6.0, 16 4 , 24 4 , 32 4 and (3 = 6.4, 32 4 , 48 4 
lattices. We study Gribov copy problem by parallel tempering (PT) gauge fixing of (3 = 2.2, 16 4 SU(2) config- 
urations. Rather unexpected features of the unique gauge (fundamental modular gauge, FMG) are observed in 
relation of physical values vs the optimizing function. 



1. The gluon propagator and the ghost 
propagator 

We measured SU(3) gluon propagator of [3 = 
6.0 on lattice 16 4 , 24 4 , 32 4 and l3 = 6.4 on 32 4 , 
48 4 [1,2]. Shown in Fig.l is the gluon dressing 
function Z(q 2 ) which is defined by the gluon prop- 
agator D lll/ (q) through 

D^(q) = (6^-^)D A (q 2 ) (1) 

as Z(q 2 ) = q 2 DA(q 2 )- Here, we adopted the logt/ 
type definition of gauge field A. Although the size 
of the sample is not large, the data scale with 
each other and consistent with the result of [3]. 
The fitted line is obtained from the model based 
on the principle of minimal sensitivity (PMS) [4,5]. 
The PMS allows one to obtain a solution that 
bridges infrared region and ultraviolet region, but 
its uniqueness is not guaranteed and in the in- 
frared region, there is discrepancy from other 
scheme that satisfies the commensurate scale rela- 
tion^]. We leave the continuation to the infrared 
region to the future study and fix the parameter y 
by the solution at the factorization scale[7,8] i.e. 
q = 1.97GeV and call it the MOM scheme. 

The ghost propagator is defined as the Fourier 
transform of an expectation value of the inverse 
Faddeev-Popov operator M. 

D$(x,y) = (tr{X a x \(M[U])- 1 \X b y)) (2) 
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Figure 1. The gluon dressing function as the 
function of the momentum q{GeV). (3 = 
6.0, 24 4 (triangle), 32 4 (diamond) and l3 = 6.4, 
48 4 (star) in logC/ version. The fitted line is 
Z(q 2 ,y) in MOM scheme. 



The lattice data scales and the PMS in MOM 
scheme fits in the region q > 0.5GeU. 



2. QCD running coupling and Kugo-Ojima 
parameter 

In the contour-improved perturbation the- 
ory [9], physical quantities 1Z are expressed in a 
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Figure 2. The ghost propagator as the function of 
the momentum q(GeV). (3 — 6.0, 24 4 (triangle), 
32 4 (diamond) and (3 — 6.4, 48 4 (star) in log t7 ver- 
sion. The fitted line is Z g (q 2 ,y) of the MOM 
scheme, which is singular at q ~ OAGeV. 
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n(q 2 )=B 1 (q 2 ) + Y / A n B n+1 (q 2 ) (3) 

71=1 

where A n is known by perturbation theory in loop 
expansion, and B n are expressed by the Lambert 
W function. The result of running coupling using 
= 237MeV fits the data of (3 = 6.4, 48 4 
data, except the peak around lGeV. 

In the j3 = 6.4, 48 4 lattice in log U version, 
we obtained the Kugo-Ojima parameter c = 
0.793(61) [10], trace e/d = 0.982(1) and the hori- 
zon function deviation h — — 0.19[11,1,2]. 

3. Gauge fixing in use of PT and Gribov 
copy problem in SU(2) 

Fundamental modular gauge (FMG: global 
maximum of the optimizing function 2^(5)) [11] is 
the unique gauge known as the legitimate Landau 
gauge without Gribov copies [12], where Fjj{g) = 
J^RetrC/^. Simulated annealing methods are 
now widely used to cope with Gribov copies in the 
maximal abelian gauge[13]. In the present work, 
we combine parallel tempering (PT) method as a 
new trial with gauge fixing for FMG, and test it in 
SU(2) gauge theory on 16 4 lattice with (3 — 2.2 



Figure 3. The running coupling a s (q) as a func- 
tion of momentum q(GeV). (3 — 6.4, 48 4 lattice. 
Fitted line is the result of the contour improved 
perturbation method. 



of the Wilson action. The PT was invented as 
one of suitable methods for simulating correctly 
low temperature equilibrium of systems possess- 
ing complex energy profile like spin glass [14, 15]. 
We do not need the exact equilibrium of Boltz- 
mann distribution Prob(g) oc exp{(3Fjj{g)} 1 but 
only the lowest energy state g of Min g E(g) effi- 
ciently where E(g) = —Fjj(g). In the PT method, 
we prepare N r replicas of site-spin (gauge trans- 
formation) system with gauge link interaction U. 
Each replica gj is given different /3, . , /3 m i n < 
f3i < Pmax- After N mc Monte Carlo steps at 
the given (3, the f3's are exchanged between adja- 
cent values according to Metropolis acceptance 
check[14,15]. During N ex iterations of the (3- 
exchange step, the lowest energy state among all 
replicas is monitored and stored. These accumu- 
lated N c gauge spins are all fixed to the Lan- 
dau gauge by the standard over-relaxation (OR) 
method (M&x x (dA) x < 10~ 4 ). We pick up the 
best copy among those, and name it as PT copy. 
In the PT method, the distribution of (3's in 
the interval [/3 m in , Anax] is technically the most 
important point so that the exchange rates are 
almost equal everywhere between replicas. We 
found the good choice for adjustme nt of the new 
/3"s from the old (3's as A/3' oc y/A/(AE/A/3) 
where A = A(3 ■ AE and A is the replica average 
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of A. Due to the actual situation of our choice, 
Pmin — 0.4 and (3 max = 30, and the system size 
of 16 4 , and N r — 24, direct application of the PT 
is practically meaningless since the Metropolis ac- 
ceptance rate is very small. We dare to shrink the 
energy scale by a factor a — 0.0002 only within 
the (3 exchange procedure to have the Metropo- 
lice acceptance rate around 0.5. Other parame- 
ters of ours are as follows, N mc = 120, N ex = 12, 
N c = 16 ~ 32. Total number of PT gauge fixed 
samples of /3 = 2.2 of SU{2) Wilson action is 67, 
and we compared them with gauge copies that 
are directly gauge fixed by OR method and are 
named as first copies. Among 67 samples, 63 
PT copies give larger Fu(l) with N c — 16 than 
first copies and the remaining 4 PT copies with 
N c = 20 ~ 32. Using these two groups of copies, 
Gribov noise is investigated in Kugo-Ojima pa- 
rameter, gluon propagator and ghost propaga- 
tor. The averaged Kugo-Ojima parameter is sup- 
pressed about 4% in PT in comparison with first 
copy. Effects of noise about 6% are seen in Fig. 4 
for the ghost propagator, but no significant ef- 
fects in the gluon propagator. For a few samples 
among 67, comparison is made for each sample 
between two groups of 16(= N c ) PT generated 
copies, and 32 copies generated(CRGT) by direct 
OR gauge fixing after random gauge transforma- 
tion, in a scatter plot of Kugo-Ojima parameter 
vs Fjj(1). It is observed in this plot that a group 
of PT copies scatter almost in a line on top-end of 
Fu(l), with all CRGT scattering in lower Ftr(l) 
area, and this fact raises questions if this "flat- 
end" feature is real one of the FMG, or if we have 
not reached the real FMG yet? This is a rather 
unexpected feature with respect to FMG. 

4. Discussion and outlook 

The momentum dependences of the gluon prop- 
agator and the ghost propagator in the infrared 
are analyzed in use of PMS method. Deviation 
of u(0) — — c from -1 (Kugo-Ojima confinement 
criterion) does not necessarily imply violation of 
colour symmetry. It could imply a question on the 
realization of BRST quartet mechanism on the 
lattice and/or even further validity of the contin- 
uum Lagrangian in the infrared region. So far as 
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Figure 4. The SU (2) colour diagonal ghost dress- 
ing function as the function of the lattice momen- 
tum squared q 2 a 2 . 1st copies(diamond) are more 
singular than PT(triangle). 



our data in off-diagonal elements of ghost prop- 
agator show, there appear no definite signals of 
colour global symmetry breaking. 

From the scatter plot the Kugo-Ojima vs 
Fu(l), it seems that our PT gauge fixing works 
well, but more careful parameter tuning and com- 
parison with simulated annealing method deserve 
a future study. Gribov noise in the ghost propa- 
gator is qualitatively similar to Cucchieri's noise 
at smaller /3 = 1.6[16]. The infrared singularity 
is weakly suppressed in FMG. 

Some revisions in numerical data of ghost and 
Kugo-Ojima parameter have been made in this 
paper due to use of conjugate gradient method. 

This work was supported by the KEK super- 
computing project No. 03-94. 
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